Internal localized eigenmodes on spin discrete breathers in antiferromagnetic chains 

with on-site easy axis anisotropy 



(N 
O 
O 



CZ2 



X3 
O 

o 



> 
m 
in 
o 
i> 
o 

(N 
O 



o 



X 



Sang Wook Kim^ and Seunghwan Kim^ 
^ Max-Planck-Institut fiir Physik komplexer Systeme, Nothnitzer Strafie 38, D-01187 Dresden, Germany 

^ Asia Pacific Center for Theoretical Physics 
Nonlinear and Complex Systems Laboratory, Department of Physics, 
Pohang University of Science and Technology, Pohang, Korea 790-784 
(February 1, 2008) 



We investigate internal localized eigenmodes of the lin- 
earized equation around spin discrete breathers in ID antifer- 
romagnets with on-site easy axis anisotropy. The threshold of 
occurrence of the internal localized eigenmodes has a typical 
structure in parameter space depending on the frequency of 
the spin discrete breather. We also performed molecular dy- 
namics simulation in order to show the validity of our linear 
analysis. 

PACS number(s): 05.45.-a, 76.50.-hg, 75.30. Ds 



Discrete breathers (DB's) are the time-periodic and 
spatially localized excitations on the nonlinear lattice 
with translational invariance which are also called 

intrinsic localized modes in the literatures The exis- 
tence of DB's was mathematically proved in the anticon- 
tinuous limit Q. Both the discreteness and the nonlin- 
earity play a crucial role in the existence of DB's. The 
spatial discreteness is quite common in nature, particu- 
larly in condensed matter physics. Recently the experi- 
mental realization of a DB was performed in anisotropic 
antiferromagnets |^ and Josephson junction ladders 

In magnetic systems, both the spin-spin exchange in- 
teraction and the on-site spin anisotropy are intrinsically 
nonlinear, so that it is quite natural to predict the ex- 
istence of DB's. Since the dissipation of spin waves in 
magnetic materials is usually weak compared with that 
of lattice vibrations in crystals, the spin lattice model has 
obvious advantages over lattice vibration models from the 
experimental point of views. Lai and Sievers have exten- 
sively studied the DB's of spin wave, namely spin DB's, 
for various situations of antiferromagnets 0. The spin 
DB's have also been recently studied in ferromagnetic 
lattices g. 

The DB can play a role of the scattering center affect- 
ing energy transport by scattering, absorbing or radiat- 
ing phonons. The scattering properties of the DB are 
closely related to the structure of the eigenmodes of the 
linearized equation around the DB itself ]9|-[T^, which 
will be called internal localized eigenmodes (ILE's). In 
particular, it was shown that the perfect transmission 
occurs at the ILE threshold, which should be tangent to 
the phonon band edge with the zero wave number 
However, when the ILE's on DB's penetrate the phonon 



band, the well known Fano resonances are obtained p2[ . 
These results can be directly applied to the spin lattice 
model considering the analogy between the lattice vibra- 
tion (phonon) and the rotation of spin (magnon or spin 
wave) . 

In this paper, we present the existence of the ILE's on 
spin DB's using the linearized equation of spin DB's in 
one-dimensional (ID) antiferromagnets with on-site easy 
axis anisotropy. The thresholds of the ILE's show some 
special structure in parameter space. The comparison 
between the prediction from linear analyses and the re- 
sults of molecular dynamics (MD) simulation will also be 
presented. 

Let us consider the antiferromagnetic chain of N spins 
with the Hamiltonian lOi 



i/-2J^S„.S„+i-i?^(5^)2 



(1) 



where the positive J and D are the exchange con- 
stant and the single ion anisotropy constant, respectively. 
Hence, the antiferromagnetic ordering and aligning with 
z direction of each spin are energetically more favorable 
in the ground state. 

Using the well known Heisenberg equation of motion 



at 



(2) 



with s,^ — {S^HS^)/S, the nonlinear equation of motion 
for s+ can be obtained in the following 



dt 
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the stationary spin DB can be obtained by using the 
following equation 
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flUJS^ = (-1)"{C( 

-I- [C(s„_i + s„+i) + - 4} 



(4) 



where T% = h/2DS and C = J/D. In this paper we 
consider only a single non- moving spin DB. The famil- 
iar dispersion relation of the extended spin wave modes. 
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TiLj = ± + 2C)2 — 4(7^ cos^ ka, can also be obtained 
by Eq. (§) by putting 
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(5) 



where a is the lattice spacing, and |/o|, |/i| <C 1. The 
linearized equation of Eq. (|^) near the spin DB for 



is given by 



ifi(n = (-l)"{C[(A„_i + - s„(B„_ii?e(e„-i) 

+ B„+ii?e(C„+i)) + An{^n-l + ^n+i) - S„(s„_i (6) 
+ s„+i)i?e(^„)] + A„^„ - s„_B„i?e(^„)} - ^W6^„, 



where An = ■\/ 1 — , i3„ = s„/2-\/l — s^, and Re{£,) is 
the real part of ^. The above Floquet (or linear) equation 
is written in the form d£,/dt = MS,. By diagonalizing 
the matrix M, the eigenmodes for small perturbation of 
spin DB's are obtained. Let us note that there exist 
two important system parameters, the relative strength 
of the exchange interaction, C{— J/D) and the frequency 
of spin DB, uji,. 
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FIG. 1. The band structure of spin waves as a function of 
C. The inset represents the amplitude of the spin oscillation 
at the central site of spin DB as a function of C for a given 
frequency of spin DB ut- From the upper they correspond to 
ujb/ujrn ~ 0.8, 0.9 1.0, 1.1 and 1.2, respectively. 

Without the exchange interaction, i.e. J = 0, the res- 
onance frequency of a local nonlinear oscillator is limited 
to LUm = since the nonlinearity of the individual sys- 
tem is soft. In the case ajb/iUm > 1, spin DB's cannot 
exist in the small coupling limit ( J — > 0), so that we ob- 
tain two distinct regions for Ub such as LOb/oJm > 1 and 
LUii/ujm < 1. In the former case (ujb/ujrn > 1) the spin 
DB's can exist only above the certain critical value of C, 
which can be explained by considering a new on-site po- 
tential modified by the neighboring sites. For simplicity 
let's consider just three spins, namely, ones on the sites 
— 1, 0, and 1, where the site corresponds to the center 
of the spin DB. In the weak coupling limit the dynamics 
of the neighboring spins at the sites ±1 can be described 



by simple harmonic oscillation with a small amplitude e, 
assuming the spin DB solution. From Eq. (^) we can ob- 
tain the following equation for the frequency of the spin 
precession at the site 



fwjSn 



(7) 



Considering e << so < 1 in the weak coupling limit, 
the above equation can be approximated to loIlo„x ~ 
2C -|- -^/l — Sq, so that the new maximum frequency 
of local nonlinear oscillators, namely w^, is given by 
LOm[2C -f 1), which means that the possible maximum 
frequency of a spin DB also increases under the coupling 
with the neighboring spins. Actually this relation exactly 
coincides with the k — 7^/2 band edge of the spin wave as 
shown in Fig. |l| in the weak coupling limit. Taking into 
account the instability caused by the resonance between 
DB's and phonons [Q, we should exclude the overlapping 
region between spin DB and the spin wave band, so that 
finally the border of the existence of spin DB's is given 
by the /c = band edge, Lu'^/Lu,n — Vl + 4C. The in- 
set of Fig. |l| shows the amplitude of the spin oscillation 
at the central site of spin DB as a function of C for a 
given frequency of spin DB Wb, where the criteria men- 
tioned above can be also confirmed numerically. In other 
words, the frequencies of spin DB's are always located 
below the spin wave band, which is nothing new in the 
sense of general criteria for the existence of a DB. We 
would like to mention, however, that this leads to some- 
what interesting consequences that a spin DB appears 
above some critical coupling strength, and disappears in 
the anticontinuous limit with the frequency of a spin DB 
fixed. The mathematical proof for the existence of DB's 
was performed in this limit jij . 
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FIG. 2. Floquet spectrum for (a) LOhjiOm ~ 0.8 (b) 
uJb/i^Om — 1.0, and (c) Ub/uJm = 1.2. See the left tick labels 
for (a) and the right tick labels for both (b) and (c). 

Figure ^ shows the evolution of the eigenvalues of the 
Floquet matrix M on unit circle (the angle (f> of the eigen- 
value e*"^) as a function of C for several w;,'s, where both 
spin wave bands and the detached branches are clearly 
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seen. For all calculations we take h = 1. The eigenmodes 
coming from the bottom of the band are observed to be 
the ones which are localized. These are ILE's mentioned 
above. The ILE's belong to either symmetric or antisym- 
metric eigenmodes with respect to their reflection sym- 
metry to the center of a spin DB. While for LUb/uJm < 1 
two ILE's appear simultaneously at C = [Fig.||(a)], for 
<^b/ijJm > 1 at first the symmetric ILE appears at a cer- 
tain value of C and then the antisymmetric one appears 
at larger C [Fig. ^(c)] . We call the former the case I and 
the latter the case II. 




FIG. 3. The threshold of occurrences of symmetric (the 
filled squares) and antisymmetric (the open squares) ILE's as 
a function of the frequency of the spin DB. Spin DB's cannot 
exist in the region under the solid line (see the text for details). 
In the inset we replot the threshold of the antisymmetric ILE's 
(the open squares) in log-log scale, which clearly shows power 
law dependence. 

Figure |^ shows the threshold of the ILE's as a func- 
tion of iOb, where the dashed vertical line Ub/ujm = 1, at 
which two spin wave bands collide each other as shown 
in Fig. 11(b), divides the parameter space into two re- 
gions for the thresholds of ILE's (i.e. the case I and II). 
The solid line represents the border of the existence of 
spin DB's described by C = [(wb/wm)^ — l]/4, below 
which the spin DB's cannot exist. It is mentioned that 
the threshold of the symmetric ILE (the filled squares) 
looks like the straight line, and that of the antisymmet- 
ric one (the open squares) shows the power law depen- 
dence, namely C oc uj'^ '^'^ as shown in the inset of Fig. ||. 
This threshold behavior was also observed in the case 
of Klein-Gordon (KG) chain with 4>'^ (double well) on- 
site potential described by V{x) = (x'' - 2x'^)/A [0, 
where LOb = 2/3, at which two neighboring phonon band 
edges with fc = collide each other, play a similar role 
as uJb/uJm = 1 in our antiferromagnet. However, this fre- 
quency does not correspond to the maximum frequency 
of the local soft nonlinear oscillator, which is ujm = 1 in 
a (j)'^ KG chain. It should also be noted that in other 
KG chains with Morse {V{x) = [1 — exp(— x)]^/2) or cu- 
bic [V{x) = — x^/3] on-site potential this is not the 



case, where for all LUb{< '-^m = 1) the cubic and the Morse 
potential correspond only to the case I and the case II, 
respectively. This has not been understood yet. 
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FIG. 4. The power spectrum of A{[^(t) with uji/curn = 1-1 
and C = 0.3. This is calculated using MD simulation with 
103 spins during 256Tb {Tt = 2TT/uJb) under the random noise 
with the intensity D — 10""" [14]. "★" indicates the symmetric 
ILE, above which the spin wave excitations are shown. 

One of the reasons why we choose this magnetic lattice 
for our investigation is that experimentally the genera- 
tion of a spin DB was reported in the quasi-lD biaxial 
antiferromagnet (C2H5NH3)2CuCl4 ^ like the system 
studied in this paper. By microwave absorption exper- 
iment the peak corresponding to the spin DB was ob- 
served in the spin wave gap. The absorption spectrum 
measured in the experiment is proportional to the imag- 
inary part of the dynamic magnetic susceptibility, which 
can be calculated using the Kubo expression j|] . 

/>oo 

A{uj)(x / dt{My{t + t')My{t'))^,e''^\ (8) 

where My{t) — X^n'^rK^) ^-^id {A)^, is the time average 
of the variable A. We calculate My{t) using MD sim- 
ulations with 103 spins starting from the solution of a 
spin DB obtained from Eq. (^ under the random ampli- 
tude noise [|l5| . Figure ^ clearly shows that an absorption 
peak, which is a several order of magnitude smaller than 
the spin DB peak, exists between the spin DB and the 
spin wave band. Without noise we can obtain only the 
spin DB peak. Although the absorption power of this 
ILE is so tiny compared with that of the spin DB, it 
should not be ignored in comparison with that of the 
spin wave band. It is worth noting that only the sym- 
metric ILE can be observed in the absorption spectrum 
since the antisymmetric ILE disappears when the sum- 
mation, X]n*ra(0> performed. Figure || shows that 
the frequency of the symmetric ILE calculated using MD 
simulation is well fitted by the prediction of the linear 
theory in the region of small C. It is remarked that in 
the experiment of Ref . ||] , somewhat broad spectrum of 
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a spin DB was observed since many spin DB's with differ- 
ent frequencies were excited simultaneously. In order to 
observe the ILE's studied in this paper, more improved 
experimental status will be needed, for example, a sin- 
gle spin DB excitation, the suppression of volume and 
surface modes, better resolution in the absorption power, 
and so on. We also note that in general the antiferromag- 
netic ordering can persist only above some critical value 
of C, roughly speaking the order of one. Even though it 
does not seem to be easy to find out the transition from 
the case I to the case II in this system in experiment 
due to this reason, it is still possible to observe the ILE's 
themselves for higher values of C. 
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FIG. 5. The frequency of the symmetric ILE given by the 
Floquet spectrum in Fig. ^ (the solid curve) and by the power 
spectrum of M^(t) using MD simulation (the filled square). 

In our previous work on discrete nonlinear Schrodinger 
(DNLS) equation (lO| we aheady noted that the DNLS 
equation is the nontrivial simplest system for studying 
scattering problems or ILE's. However, it has only a 
single parameter since the change of the breather fre- 
quency ujh can be compensated by scaling SI — > Q./uJi, 
and C C/oJh, where Q. and C are the frequency of 
the small perturbation and the linear coupling strength, 
respectively. We would like to point out that the ID an- 
tiferromagnet with easy axis anisotropy provides a good 
intermediate example with two spin wave bands (or scat- 
tering channels) and two parameters (C and u>h) between 
the case of DNLS equation with two phonon bands and 
one parameter (C/wb) and that of the KG chain with an 
infinite number of phonon bands and two parameters. 

In summary we have investigated the characteristics 
and the possibility of experimental observation of the 
ILE's in ID antiferromagnets with easy axis anisotropy 
using both linear analyses and MD simulations. The 
structure of the thresholds of the ILE's strongly depends 
on the frequency of a spin DB. For Wb/wm < 1 the ILE's 
appear at C = simultaneously, while for LOb/uim > 1 the 
symmetric ILE occurs at first at the certain value of C 
and then the antisymmetric one appears at the larger C. 



It is shown that the frequency of the ILE's calculated us- 
ing MD simulation is well fitted by the prediction of the 
linear theory. We hope that these ILE's will be observed 
in experiments. 
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thanks Mikhail Fistul for helpful discussions. 
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